Application Serial No. 10/670,1 16 

Amendment Dated August 16, 2004 

Reply to Notice of Allowance of May 14, 2004 

AMENDMENTS TO THE SPECIFICATION 
Please replace the paragraph at page 1, lines 10-15 with the following rewritten paragraph: 
The present invention relates to estimating the precise position of a stationary or moving 
object using multiple satellite signals and a network of multiple receivers. The present invention 
is particularly suited to position estimation in real-time kinetic environments where it is desirable 
to take into account the spatial distribution of the ionosphere delay. 

Please replace the paragraph at page 1, lines 17-28 with the following rewritten paragraph: 

Satellite navigation systems, such as GPS (USA) and GLONASS (Russia), are intended for 
accuracy accurate self-positioning of different users possessing special navigation receivers. A 
navigation receiver receives and processes radio signals broadcast by satellites located within 
line-of-sight distance, and from this, computes the position of the receiver within a pre defined 
predefined coordinate system. However, for military reasons, the most accurate parts of these 
satellite signals are encrypted with codes only known to military users. Civilian users cannot 
access the most accurate parts of the satellite signals, which makes it difficult for civilian users to 
achieve accurate results. In addition, there are sources of noise and error that degrade the 
accuracy of the satellite signals, and consequently reduce the accuracy of computed values of 
position. Such sources include carrier ambiguities, receiver time offsets, and atmospheric effects 
on the satellite signals. 

Please replace the paragraph at page 2, lines 10-27 with the following rewritten paragraph: 

In a first aspect of the present invention, an exemplary apparatus/method comprises 
receiving the known locations of a first base station and a second base station, obtaining a time 
offset representative of the time difference between the clocks of the first and second base 
stations, and obtaining measured satellite data as received by the rover, the first base station, and 
the second base station. The measured satellite data comprises pseudo-range information. The 
exemplary apparatus/method generates a first set of residuals of differential navigation equations 
associated with a set of measured pseudo-ranges related to a first baseline (R-Bl) between the 
rover and the first base station. The residuals are related to the measured satellite data received 
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by the rover station and the first base station, the locations of the satellites, and the locations of 
the rover station and the first base station. The exemplary apparatus/method also generates a 
second set of residuals of differential navigation equations associated with a set of measure 
pseudo-ranges related to a second baseline (R-B2) between the rover and the second base 
station. These residuals are related to the measured satellite data received by the rover station 
and the second base station, the locations of the satellites, and the locations of the rover station 
and the second base station. The exemplary apparatus/method estimates the rover's location 
from the first set of residuals, the second set of residuals, and the time offset between the clocks 
of the first and second base stations. 

Please replace the paragraph at page 3, lines 7-1 7 with the following rewritten paragraph: 

In a fourth aspect of the present invention, which may be applied with any of the aspects 
of the present invention described herein, the exemplary apparatus/method obtains measured 
satellite carrier phase data as received by the rover and the first base station (the first base lin e 
baseline) . The exemplary apparatus/method generates a fourth set of residuals of differential 
navigation equations associated with a set of measured carrier phase data related to the first base 
Hne baseline , and resolves the cycle ambiguities in the carrier phase data from the fourth 
residual and one or more of the first, second, and third residuals. The resolved cycle ambiguities 
may take the form of floating ambiguities, fixed-integer ambiguities, and/or integer ambiguities. 
The exemplary apparatus/method estimates the rover's location further with the fourth set of 
residuals and the resolved cycle ambiguities associated with the first base lin e baseline . 

Please replace the paragraph at page 3, lines 18-32 with the following rewritten paragraph: 

In a fifth aspect of the present invention, which may be applied with any of the aspects of 
the present invention described herein, the exemplary apparatus/method obtains measured 
satellite carrier phase data as received by the second base station, and further obtains a set of 
cycle ambiguities related to a set of satellite phase measurements associated with the baseline 
between the first and second base stations (B1-B2). The exemplary apparatus/method generates 
a fifth set of residuals of differential navigation equations associated with a set of measured 
carrier phase data related to the second base line between the rover and the second base station 
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(R-B2), and resolves the cycle ambiguities in the carrier phase data from the fifth residual, the set 
of cycle ambiguities associated with the baseline between the first and second base stations, and 
one or more of the first, second, third, and fourth residuals. The resolved cycle ambiguities may 
take the form of floating ambiguities, fixed-integer ambiguities, and/or integer ambiguities. The 
exemplary apparatus/method estimates the rover's location further with the fifth set of residuals 
and the resolved cycle ambiguities associated with the second base lin e baseline (R-B2). 

Please replace the paragraph at page 4, lines 15-27 with the following rewritten paragraph: 

In a seventh aspect of the present invention, which may be applied with any of the 
aspects of the present invention described herein, the exemplary apparatus/method obtains a first 
set of first ionosphere delay differentials associated with the satellite signals received along the 
bas e line baseline formed by the first and second base stations, and generates ionosphere delay 
corrections to one or more of the above-described first through sixth residuals from the first set 
of first ionosphere delay differentials, the locations of the first and second base stations, and an 
estimated location of the rover station. As a further option, the exemplary apparatus/method 
obtains a set of second ionosphere delay differentials associated with the satellite signals 
received along the base line baseline formed by the first and third base stations (or the base line 
baseline formed by the second and third base stations), and generates the ionosphere delay 
corrections to the one or more of the above-described first through sixth residuals further from 
the second set of ionosphere delay differentials and the location of the third base station. 

Please replace the paragraph at page 4, line 28-32 with the following rewritten paragraph: 

In an [[eight]] eighth aspect of the present invention, which may be applied to the 
seventh and further aspects of the present invention described herein, the exemplary 
apparatus/method forms one or more of the first through sixth residuals to account for second 
order effects in the ionosphere delay corrections applied to the baselines associated with the 
rover, and further generates an estimate of the second order effects. 
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Please replace the paragraph at page 5, lines 1 7-24 with the following rewritten paragraph: 

In an eleventh aspect of the present invention, which may be applied with any aspects of 
the present invention described herein which account for ionosphere delays, the exemplary 
apparatus/method generates the first set of first ionosphere delay differentials, and optionally the 
second set of first ionosphere delay differentials. As a further option, the apparatus/method 
generates a third set of ionosphere delay differentials associated with the satellite signals 
received along the base lin e baseline formed by the second and third base stations, and generates 
therefrom three sets of ionosphere delay differentials which are s elf consi s tent self-consistent . 

Please replace the paragraph starting at page 6, line 19 and ending at page 7, line 4 with the 
following rewritten paragraph: 

FIG. 1 is a perspective view of a rover station (R) and three base stations (Bl, B2, B3) in 
an exemplary network according exemplary embodiments of the present invention, and FIG. 2 is 
a top-plan schematic drawing thereof The present invention pertains to estimating the position 
of the rover station with information provided from two or three of the base stations and with 
satellite measurements made by the rover station. Referring to FIG. 1, each station has a receiver 
that receives the satellite positioning signals with a satellite antenna (shown as a substantially flat 
disk). A plurality of satellites S1-S4 are depicted in FIG.l, with the range between each satellite 
and each antenna being depicted by a respective dashed line. In the example shown [[and]] in 
FIG. 1, the Rover station is operated by a human user, and has a positioning pole for positioning 
the Rover's satellite antenna over a location whose coordinates are to be determined. The 
Rover's satellite antenna is coupled to the receiver's processor, which may be disposed in [[the]] 
a backpack and carried by the user. The user interacts with the receiver's processor through a 
keypad/display. The receiver also has a radio modem (more formally a demodulator) that can 
receive data from the base stations through a conventional RF antenna and relay it to the 
processor. In other embodiments, the Rover's satellite antenna may be mounted to a vehicle, and 
the receiver may operate independently and automatically without the need of a human operator. 
The present invention is applicable to these and other physical embodiments. 
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Please replace the paragraph at page 7, lines 5-16 with the following rewritten paragraph: 

The satellite signals comprise carrier signals which are modulated by pseudo-random 
binary codes, which are then used to measure the delay relative to a local reference clock or 
oscillator. These measurements enable one to determine the so-called pseudo-ranges between 
the receiver and the satellites. The pseudo-ranges are different from true ranges (distances) 
between the receiver and the satellites due to variations in the time scales of the satellites and 
receiver and various noise sources. To produce these time scales, each satellite has its own on- 
board atomic clock, and the receiver has its own on-board clock, which usually comprises a 
quartz crystal. If the number of satellites is large enough (four or more), then the measured 
pseudo-ranges can be processed to determine the user location (e.g., X, Y, and Z coordinates) 
and to reconcile the variations in the time scales. Finding the user location by this process is 
often referred to as solving a navigational problem or task. 

Please replace the paragraph at page 7, lines 17-27 with the following rewritten paragraph: 
More specifically, the GPS system employs a constellation of satellites in orbit around the 
earth at an altitude of approximately 26,000 km. Each GPS satellite transmits microwave radio 
signals in two frequency bands located around 1575.42 MHz and 1227.6 MHz, referred to as LI 
band and L2 band, respectively. The GPS LI -band signal is modulated by a coarse/acquisition 
code (C/A) and a precision ranging code (P-code). The L2-band signal is binary modulated by the 
P-code. The GPS C/A code is a pseudo-random (PR) Gold code that is specific to each satellite and 
is used to identify the source of a received signal. The P-code is a pseudo-random code signals 
signal and is also specific to each satellite, having a symbol rate which is ten time more than C/A, 
which reduces the granularity by a factor of ten. The GPS satellite transmission standards are set in 
detail by the ICD-GPS-200, Revision Q ARINC Research Corporation, 10 October, 1993. 

Please replace the paragraph at page 8, lines 5-23 with the following rewritten paragraph: 
The distance between a receiver and a satellite (called the "receiver-to-satellite range") is 
determined by measuring the time that it takes for the signal to pass from the satellite to the receiver, 
provided that the position of the satellite is known. The satellites and receivers have internal clocks 
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that are synchronized to a single GPS time. For each satellite signal being tracked, the receiver 
generates a local version of the satellite's expected PR-code, and then retards that version in time 
until the local version correlates (Le., matches) with the received satellite signal. Thereafter, the 
satellite signal is tracked by advancing or retarding the local version of the PR-code. The carrier 
phase of the satellite can also be tracked, which is usually done by tracking the Doppler shift of the 
satellite signal. The positions of the satellites are, except for minor variations, highly predictable as 
a function of time, and the receiver generally carriers a model of the satellite's position as a function 
of GPS time. In theory, by determining the ranges to three different satellites, the receiver can 
perform a three-dimensional triangulation to find its position. But because of limitations in the 
accuracy of the receiver's clock, the internally generated time is offset somewhat from true GPS 
time. Thus, the ranges to at least four different satellites are simultaneously measured in order to be 
able to solve for four unknowns, namely the three coordinates of the position of the receiver 
location (e.g., x, y, and z) and an offset of the receiver clock time from the GPS time. The location 
is usually performed with respect to the defined Cartesian coordinates frame. 

Please replace the paragraph starting at page 8, line 30 and ending at page 9 y line 9 with the 
following rewritten paragraph: 

The desire to guarantee the solution of navigational tasks with accuracy better than 10 
meters, and the desire to raise the stability and reliability of measurements, have led to the 
development of the mode of "differential navigation ranging," also called "differential 
navigation" (DN). In the DN mode, the task of finding the user position is performed relative to 
a Base station (Base), the coordinates of which are known with [[the]] high accuracy and 
precision. The Base station has a navigation receiver that receives the signals of the satellites 
and processes them to generate measurements. The results of these measurements enable one to 
calculate corrections, which are then transmitted to a roving GPS receiver, which the user has set 
up. We call this GPS receiver the "Rover station," or "Rover receiver." By using these 
corrections, the roving GPS receiver gains the ability to compensate for the major part of the 
strongly correlated errors in the measured pseudo-ranges, and to substantially improve the 
accuracy of the estimate of its position. 
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Please replace the paragraph at page 12, lines 13-21 with the following rewritten paragraph: 

Even with all of th e abov e described the above-described processing to account for 
clock offsets and carrier ambiguities, there are additional factors which affect the accuracy of 
measurements made with GPS and/or GLONASS signals. As one factor, the trajectory of each 
satellite (or its initial data called the ephemeris), is elliptical and is affected by natural causes such as 
solar winds. The accuracy of any measurement is dependent upon knowledge of the position of the 
satellites at a certain time. An estimate of the ephemeris is calculated on earth for each of the 
satellites and is periodically uploaded to the satellite. The position information of a satellite is 
encoded onto a low frequency (50 Hz) signal which is modulated on to one of the carrier signals, 
and transmitted to the GPS receiver on earth. 

Please replace the paragraph at page 12, lines 22-33 with the following rewritten paragraph: 
Two additional and important factors that affect the accuracy of measurements made with 
GPS and/or GLONASS signals are the effects of the troposphere and ionosphere on propagation of 
signals from the satellites to the receivers. The troposphere is the lower part of the atmosphere and 
variations in the temperature, pressure, and humidity lead to spatial variations in the signal 
propagation. The ionosphere is at the upper part of the atmosphere and has a slice of ionized gas at 
the altitude of around 300 km. The density of ionized particles is sufficiently high to affect the 
propagation of electromagnetic signals, and has a spatial variation and a time variation. The 
ionosphere effect becomes even more important during years of high solar activity. These 
variations in the troposphere and ionosphere introduce variable delays in the propagation of the 
satellite signals to the receivers. If the base and rover are widely separated, these delays will be 
significantly different for the base and rover stations, and will introduce error into the estimation of 
the rover's position. 

Please replace the paragraph at page 13, lines 1-12 with the following rewritten paragraph: 

With that general overview, we now describe the invention in greater detail. For the sake 
of simplifying the presentation, and without loss of generality, we assume "AT satellites 
identified by the index u s" s = 1, . . N, and we assume that each receiver can track the satellite's 
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LI band signal and the s atellite s satellite's L2 band signal. The receivers extract timing 
information from the satellite signals, and report this information [[a]] as predefined increments 
k of time, which we call epochs k. The time between epochs can be selected by the user, and 
generally ranges between 0.1 seconds to 2 seconds, with 1 second being typical. The clocks of 
all of the receivers are typically accurate to within several milliseconds of the true GPS time, and 
for practical purposes the receivers can determin e d determine the number k of the current 
epoch from their clocks. The following timing information can be extracted from each satellite 
"s" at each epoch *T' by each receiver V: 

Please replace the paragraph starting at page 16, line 10 and ending at page 17 y line 2 with the 
following rewritten paragraph: 

The difference form of [3 A] is generated by forming two in s tanc e instances of form [1A] for 
receivers V and "q", and then subtracting the two instances (the instance for receiver "q" is 
s ubtract subtracted from the instance for receiver "r"). Difference forms [3B], [3C], and [3D], 
are formed in a similar manner from corresponding instances of forms [IB], [1C], and [ID], 
respectively. Forms [3A]-[3D] are the single-differences of the navigation equations. The 
benefit of forming the between-station single differences is that the error term representing the 
time offset of the satellite clock, r/, is cancelled out in the differences. We emphasize that the 
forms [3A]-[3D] can be applied to each satellite that can be observed by receivers "q" and "r." 
Higher-order differences of the navigation equations, such as the double-differences of the 
navigation equations, can be formed and are known to the art. For example, a common double- 
difference equation is the difference between two s ingle s- diff e r e nce single-difference equations 
associated with a common baseline, but with each single-difference equation being based on a 
different satellite. The present invention may also use these higher-order differences of the 
navigation equations, although the single differences of the navigation equations are currently 
preferred. Each of these difference forms is generically referred to as a differential navigation 
equation. 
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Please replace the paragraph at page 17, lines 5-26 with the following rewritten paragraph: 
Because the locations of the base stations are known, and because the locations of the 
satellites as a function of epoch time (k) are also known, the between station between-station 

differences of each true range A qr R s k = Kk~ R qk associated with each satellite 'V can be 

generated from existing information. If vectors X r , X q , and X s represent the locations of receiver 
V, receiver "q'\ and satellite V, respectively, in a coordinate system (e.g., a Cartesian 
coordinate system), and the operator || * || repre s ent represents the distance operator in that 

coordinate system, then A qr R s k = II X s - x r II - II X s - Xq II- This reduces by one the number of 

unknowns in forms [3A]-[3D]. In addition, the Goad-Goodman model may be used to model the 
difference in troposph e r e tropospheric effects (with an error less than a few percent), and thus 

the difference A^ r 7^ may be estimated based on the positions of the receivers and the satellite 

(more specifically, the angles between the receivers and s atellite s satellite) . The noise sources 

A q r n LliS , A qr n L2 ' s , A q r v Ll ' s ,-and A qr v L2,s can never be known, but they are generally 

zero-mean and their effect can be reduced by averaging. Thus, the number of solvable 

unknowns in forms [3A]-[3D] may be reduced to the following six: (1) A qr r k , (2) A qr I s k , 

(3) A qjr N Ll > s 9 (4) A q9r i// Ll , (5) A q r N L2 > s . and (6) A^y/ 12 . The first unknown varies 

with time and is common to all of the satellites being tracked by receivers u q" and "r The 
second unknown varies with time and is specific to the satellite V being tracked by the pair of 
receivers "q" and "r." The third and fifth unknowns are each specific to the satellite V being 
tracked by the pair of receivers and V, and each does not normally vary with time unless a 
cycle slip occurs in the receiv e rs receiver^ phase-lock loop. The [[forth]] fourth and sixth 
unknowns are specific to the pair of receivers "q" and "r," and do not normally vary with time. 

Please replace the paragraph starting at page 18, line 27 and ending at page 19, line 4 with the 
following rewritten paragraph: 

Forms [4A]-[4D] comprise the known terms of the single-differences of the navigation 
equations, and we call them the "residuals" of the single differences, or more generally the 
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"residuals of a set of differential navigation equations/' The magnitude of each residual is 
generally less than the magnitudes of one or more of the known terms that [[forms]] form the 
residual. 

Please replace the paragraph at page 19, lines 5-70 with the following rewritten paragraph: 

Then, the unknowns can be estimated from the residuals using the following forms by 
any of a number of solution methods known to the art: 



A q,r b k 2 ' $ = ~ cA q,r T k + 



Please replace the paragraph at page 22, lines 5-8 with the following rewritten paragraph: 

The solvable unknowns for form [9] includ e s include 8X k , as well as A^q*"* > ^qfi J k > 

Ag i0 N LlyS , and A q 0 N L2,s . As described below in greater detail, estimation processes for these 

unknowns generally include instances of forms [9A]-[9D] for several satellites, and for several 
epochs of time. 

Please replace the paragraph at page 24, lines 7-20 with the following rewritten paragraph: 

According to general aspects of the present invention, the position of the Rover is 
estimated by forming a primary baseline between the Rover and one of the base stations, usually 
the closest base station, and then forming one or more secondary baselines from the Rover to 
other base stations. Then, one or more of the above relationships are applied to the configuration 
of stations to relate the measured data associated with the secondary baseline(s) with the 
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measured data associated with the primary baselin e baseline, using one or more of the 
unknowns associated with the baselines between base stations through forms [11]. To simplify 
the presentation of the present invention, and without loss of generality, we will use base station 
Bl to form the primary baseline with the Rover, and we use base stations B2 and B3 to form the 
secondary baselines with the Rover. Varying amounts of measured data from the secondary 
baselines may be related with the primary baseline. In general, the greater amount of 
information so related increases the accuracy of the estimated position of the rover, and/or 
enables a greater spacing between base stations. 

Please replace the paragraph at page 25, lines 10-19 with the following rewritten paragraph: 

The terms in forms [12A], [14A], and [16A] are the residuals of differential navigation equations 
and are known. The unknowns are contained in the right-hand sides of forms [13A], [15 A], and 
[17A]. There are a total of six solvable unknowns: A^qT^ > &2$ T k > &3,0 T k > m & ^ e three 

components of 5X k . However, using form [11 A], A 2) o r £ may be related to Aq^ through the 
estimated base-station time offset A2 5 i^ as follows: &2$ r k = ^2 9 l r k^^0,l T k = ^2 y \ T k + ^\,0 T k • 
In a similar manner, ^fi T k may be related to Aq ^ through the estimated base-station time 
offset ^\ z k as follows: A^qZ^ = ^^k^^Q^k^ A3 )1 r^ + A 1) o^ • These two forms reduce 
the number of true solvable unknowns to four, and the following modified set of forms may be 
used to estimate the true solvable unknowns: » 

Please replace the paragraph at page 25, lines 25-27 with the following rewritten paragraph: 

As an equivalent, one may view forms &2,0 r k = ^2,\ T k + ^\,0 T k m ^ ^3fi T k = ^3 9 l T k + ^l 9 0 T k 
as increasing the total number of forms by two, and may then estimate the six unknowns of the 
expanded form set: 
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Please replace the paragraph at page 26, lines 4-8 with the following rewritten paragraph: 

The rover's location (X 0 = X 0 + SX k ) may be estimated from each of the above form 

sets using any one of several [[know]] known methods. Here, we demonstrate estimating the 
rover's location with a least squares fitting approach based on the first form set. After generating 

the residuals A^j 1,5 f° r N satellites (s=l to s=N) for the baseline between the rover and the 
first base station (R-Bl), we collect them together into a first vector: 

Please replace the paragraph starting at page 27, line 13 and ending at page 28, line 2 with the 
following rewritten paragraph: 

For the reader not familiar with the GPS art, we briefly note that matrix C n k generally 

comprises a diagonal matrix, with each diagonal element being related to the noise sources in the 
two receivers that define the baseline (in this case the rover and one of the base stations). A 
covariance factor is usually associated with each satellite signal received by each receiver, and 
this covariance factor is usually related to the signal-to-noise ratio of the signal (as received by 
the receiver) and the elevation angle of the satellite (the multipath error has a strong correlation 

with the elevation angle). Each diagonal entry of matrix C n k usually comprises an addition of 

the two covariance factors associated with the two receivers which contributed contribute to 

M s 

the underlying noise quantity, e.g., base station Bl and the rover for noise quantity A 10 "^ ' . 
For more details and information on the generation of the covariance matrices, the reader is 
direct directed to A. Leick, GPS Satellite Surveying, John Wiley & Sons, (1995). Matrix A* 
is known as the observation matrix, and it relates the solvable unknowns to the residuals. The 
solvable unknowns may also be estimated by other processes, such as various Kalman filtering 
processes. 
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Please replace the paragraph starting at page 28, line 15 and ending at page 29, line 5 with the 
following rewritten paragraph: 

Thus, the above exemplary embodiments provide methods of estimating the location of 
the rover station (R) with the use of a first base station (Bl) a second base station (B2), and 
optionally a third base station (B3) or more base stations. In summary, the locations of the base 
stations were obtained, and one of the base stations (e.g., Bl) was selected to form a primary 
baseline with the Rover station. We refer to this base station as the primary base station, and the 
other base stations as the secondary base stations. Additionally, for each of the secondary base 
stations, the time offset representative of the time difference between the clocks of the primary 
station and secondary base station is obtained. Also, measured satellite data as received by the 
rover, the primary base station, and the secondary base station(s) is obtained. From this, the set 

of residuals A^B^ 1 associated with the primary baseline (R-Bl) is generated, and is related to 

the measured satellite data received by the rover station and the first base station, the locations of 
the satellites, and the locations of the rover station and the first base station. Similarly, the set(s) 

of residuals A 2j oB^ 1 , A 3 qB^ 1 , etc. associated with the secondary baseline(s) are generated, 

each set of residuals being related to the measured satellite data received by the rover station and 
the secondary base station, the locations of the satellites, and the locations of the rover station 
and the secondary base station. Thereafter, the rover's location is estimated from the above sets 
of residuals, the time offset between the clocks of the base stations, and typically an observation 
matrix. 

Please replace the paragraph at page 31, line 14 and ending at page 32, line 2 with the 
following rewritten paragraph: 

The terms in forms [12C, D], [14C, D], and [16C, D] are the residuals of the differential 
navigation equations and are known. The unknowns are contained on the right-hand sides of 
forms [13C, D], [15C, D], and [17C, D]. For N satellites, there are 6*N equations and a total of 
(6+ 6*jV) solvable unknowns: A 10 r^ , A 2) o*> , A 3 0 ^ > the three components of 5X k , and N 

instances of each of A\ 9 oN Lls , A 2 ,0^ Lls > ^3fi^ LU > A \fi^ L2,S > A 2,0& L2 ' S > ^ d 
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A 30 N L2)S for s = 1 to s = N satellites. However, using forms [1 1 A], [1 1C], and [1 ID], the 
number of unknowns may be reduced to 4+2 *N. As we saw above, A 2 fiT k may be related to 
&0,\ T k through the estimated base-station time offset A 2j ir^ using form [1 1 A] as follows: 
&2,0 T k = A 2,l r A:"" A 0,l r A: = A 2,l r A: + A l,0 r A: • In a similar manner, A 3 0 ^ [[wasll may be related 
to A 0j ir£ through the estimated base-station time offset ^\ T k as follows: A^qt^ = A^t^- 
A 0 iT k = A 31 r^ + Ai o^jfc • Form t 1 1C] may be used to related relate each of the LI band 
ambiguities of the secondary baselines to the primary baseline as follows: 

Please replace the paragraph at page 34, lines 11-15 with the following rewritten paragraph: 

The solvable unknowns [<HT t , A 10 r fc , A 10 N L1 , A 10 N L2 ] may be estimated by applying a 

number of processes to the combination of forms [21] and [18+]. Further below, in a subsequent 
section, we described describe a preferred process that can be used on the combination of forms 
[21] and [18+]. Here, we describe how a least squares process may be applied to combined 
forms of [18+] and [21]. 

Please replace the paragraph at page 35, lines 2-4 with the following rewritten paragraph: 
A least squares process is applied to form [23] to generate the floating ambiguities over several 
epochs, and the floating ambiguities are averaged to generate [[a]] an estimate of the floating 
ambiguities. For each epoch, the following may be generated: 

Please replace the paragraph at page 37 lines 7-11 with the following rewritten paragraph: 

where A 24* is the observation matrix of form [24] (the matrix which is multiplying the 

solvable unknowns in form [24]). In form [25], the floating ambiguities may be used in place of 
the fixed-integer ambiguities. However, lower accuracy generally results, although the 
estimation speed is increased since the step of generating the fix - integ e r fixed-integer 
ambiguities from the floating ambiguities may be omitted. 
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Please replace the paragraph starting at page 37, line 12 and ending at page 38, line 8 with the 
following rewritten paragraph: 

With the large amount of measurement data afforded by the present invention, various 
residuals may be omitted from use in the estimation process. For example, we may [[only]] 
work only with LI -band data of all three stations and [[only]] use only the residuals associated 
with this data. We may also [[only]] work only with data and residuals from two base stations 
(the primary baseline and one secondary baseline). For applications needing lower accuracy, we 
may work with the phase and ps e udorang e pseudo-range data of the primary baseline (R-Bl) 
and just the p s eudorange pseudo-range data or phase data of one of the secondary baselines. 
Also, one can undertake an analysis of the satellite constellation and select the satellites which 
should provide the highest accuracy, using the pseudo-range and phase data from the LI- and 
L2- bands. Finally, while we have illustrated the invention with the single-difference navigation 
equations, it may be appreciated that higher-order differences of the navigation equations may be 
used. Such higher-order differences have known quantities (which form the residuals) and 
similar unknowns to be solved for, which can be solved by the above-described methods. 

Please replace the paragraph at page 39, lines 2-20 with the following rewritten paragraph: 

In the above embodiments, the ionosphere delays were assumed to equally affect the 
rover and base stations, and the between-station differences were neglected. When the base 
stations and rover are separated by large distances, better accuracy may be obtained by taking 
into account the ionosphere delays. This can be done in a number of ways according to the 
present invention. FIG. 3 provides a representation of the ionosphere delays of one satellite "s" 
at the base and rover stations. Shown is a 3-d Cartesian system having two planar axes, north (n) 
and east (e), to represent the terrain on which the stations are located, and a vertical axis to 
represent the ionosphere delays of a satellite "s" as a function of the terrain. The ionosphere 
delay will be different for each satellite. The locations of the rover station R and three base 
stations Bl, B2, and B3 are indicated in the north-east plane of the figure. The ionosphere delays 

for each of these stations are indicated as Ifa , ifa , m & H t k > respectively. The preferred 
embodiments of the processing of the base station data, which is more fully described below, 
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generates estimates for the between base s tation s between-base station differences in the 
ionosphere delays A^l, Ai3/|,and A23/I , which we call ionosphere delay differentials. 
From two of these differentials, an estimate for the ionosphere delay between the rover and any 
of the base stations can be generated. Here, we are most interested in the difference A 1? q^ 
associated with the primary baseline, the one between the rover and the first base station. 

Please replace the paragraph starting at page 43, line 18 and ending at page 44, line 22 with 
the following rewritten paragraph: 

Having thus described three general groups of embodiments, we now describe an 
exemplary rover station 100 in FIG. 5 that may be used to implement any of the above-described 
embodiments. De s cription A description of rover station 100 is provided in conjunction with a 
flow diagram shown in FIG. 6. Referring to FIG. 5, rover 100 comprises a GPS antenna 101 for 
receiving navigation satellite signals, [[and]] an RF antenna 102 for receiving information from 
the base stations, a main processor 1 10, an instruction memory 112 [[and]] a data memory 114 
for processor 110, and a keyboard/display 1 15 for interfacing with a human user. Memories 112 
and 1 14 may be separate, or diff e r e nc e different sections of the same memory bank. Rover 100 

further comprises a satellite-signal demodulator 120 for generating the navigation data p^ s , 

P§\ S > $0 Y ' m ^ ^0 A S ^ or ea °k e P 0C h * fr° m * e signals received by GPS antenna 101, which 
is provided to processor 110. Rover 100 also comprises a base-station information demodulator 
130 that receives information signals from the base stations by way of RF antenna 102. 
Demodulators 120 and 130 may be of any conventional design. The information received by 
demodulator 130 includes the positions (Xi, X 2 , X 3 ) of the base stations, the satellite- 
navigational data (e.g., k, p^\ s , p^ s , <f>fy s , and , r =1,2,3) received by each base 

station at each epoch and the between-base station unknowns (e.g., A r q r k , A r q N u,s , A r q y/ U , 

etc., {r,q} = {1,2}, {2,3}, {3,1}). Each set of information may be transmitted on a respective 
frequency channel. The between-base station unknowns may be generated by the first base 
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station Bl, and thereafter transmitted to the rover by the first base station. The first base station 
may receive the satellite-navigation data from the other base stations in order to compute the 
between-base station unknowns. Methods of generating the between-base station unknowns are 
described below in greater detail in the section entitled Between Base Station Processing, Part IL 
As another approach, the betw e en th e between-base station unknowns may be generated by 
rover 100 locally by a between-base station processor 140, which receives the positions of the 
base stations and their satellite navigation data from base station information to modulator 130. 
Processor 140 may implement the same methods described in greater detail in the below section 
entitled Between Base Station Processing, Part II . In addition, processor 140 may comprise its 
own instruction and data memory, or may be implemented as part of main processor 110, such as 
by being implemented as a sub-process executed by main processor 110. 

Please replace the paragraph starting at page 44, line 23 and ending at page 45, line 5 with the 
following rewritten paragraph: 

Main processor 110 may be configured to implement any and abov e de s cribed of the 
above-described embodiments by the instructions stored in instruction memory 112. We 
describe implementation of these embodiments with respect to FIG. 6, where certain of the steps 
may be omitted when not needed by a particular embodiment. In step 202, the locations Xi, X 2 , 
X 3 of the base stations Bl, B2, B3 are received by base-station information modulator 130 and 
conveyed to main processor 1 10. These locations, and the location of the rover station, are 
measured at the phase centers of the GPS antennas. Thus, RF antenna 102 and demodulator 130 
provid es provide means for receiving the locations of the first base station and the second base 
station. Also in step 202, main processor 110 generates an initial estimated location for rover 
100, which is relatively course coarse . This may be the center of the triangle formed by the base 
stations, or may be derived from a conventional s ingle point single-point GPS measurements (as 
opposed to a differential GPS measurement), or [[maybe]] may be generated by other means. 
The means for performing this step is provided by main processor 1 10 under the direction of an 
instruction set stored in memory 112. 
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Please replace the paragraph at page 45, lines 14-27 with the following rewritten paragraph: 
In step 206, main processor 1 10 is provided with the measured satellite-navigation data 
(e.g., pj:Y , , , and <j)^ s ) as received by the rover (r=0), the primary base 

station(r=l), and the secondary base station(s) (r=2,3) at one or more time moments k (epochs). 
The data is provided by the modulators 120 (for the Rover data) and 130 (for the base station 
data). Although these sets of data for each epoch k may be received at slightly different times 
(because the base stations are at different distances from the Rover), the data sets are time- 
stamped with the epoch identifier (which is conventional practice), and can be stored in a 
synchronized queue until all the data sets for the epoch are received. Also during this step, main 
processor 110 determines the positions of the satellites for these time moments from orbital 
predictions, and generate generates the computed ranges of the satellites to the rover and base 
stations based on the positions of the satellites and the stations. The means for performing this . 
step is provided by main processor 110 under the direction of an instruction set stored in memory 
1 12, with the computed information being stored in data memory 114. 

Please replace the paragraph at page 46, lines 5-18 with the following rewritten paragraph; 

The next step 210 is optional, depending upon the embodiment being implemented. In 
this step, rover 100 obtains, for each secondary base station, a set of satellite-phase cycle 
ambiguities related to the baseline between the secondary and primary base stations for one or 
more of the frequency bands (e.g., LI and L2) at the time moments k. These ambiguities may be 

in floating form (e.g, A^A^ 1 ' 5 , A 2i \N L2i \ A 3fl N Ll > s , A 31 iV Z ' 2 ' 5 ), fixed-integer form (e.g, 
A 2,\N LliS > A 2,\N L2 ' S i A 3,1^ Z1 ' 5 > A 3,l^ 12 ' 5 ) or integer plus fractional phase form (e.g., 
A 2 ^, A 2jl ^A 2?1 A^, A 2ilV , L2 ; A 3yl N Ll > s , A^ 1 ; A 3jl tf 12 '*, A 3jl ^ 2 ). 

Rover 100 may obtain this information by directly receiving it from the primary base station Bl 
by way of demodulator 130, or may obtain [[at]] this information by using between-base station 
processor 140 to generate it, as indicated above. Either of these approaches provides the means 
for obtaining this information. As used herein, the term "obtain" encompasses both the receiving 
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of the information from an outside source (e.g., the primary base station) and the generating the 
information locally with processor 140. 

Please replace the paragraph at page 47, lines 4-7 with the following rewritten paragraph: 

As an optional step 216, main processor 110 generates ionosphere corrections and adds 
these corrections to the residuals. The means for performing this step is provided by main 
processor 110 under the direction of an instruction set stored in memory 112, with the 
corrections being stored in data memory 1 14. 

Please replace the paragraph at page 47, lines 8-15 with the following rewritten paragraph: 

Finally, in step 218, main processor 110 estimates e s timate the rover's location at the one 
or more time moments k from the sets of residuals, the time offset between the clocks of the base 
stations, the sets of satellite-phase cycle ambiguities related to the baseline between the primary 
base station and the secondary base stations (optional), and an observation matrix. The 
estimation may be done by the previously-described methods. The means for performing this 
estimation step is provided by main processor 110 under the direction of an instruction set stored 
in memory 1 12, with the computed information being stored in data memory 1 14. 

Please replace the paragraph at page 47, lines 16-20 with the following rewritten paragraph: 

Keypad/display 115 may be used to receive an instruction from the human user to 
commence an estimation of the position of the Rover, and to provide an indication of the 
estimated position of the Rover to the user. For some applications, it may be appreciated that 
human [[and]] interaction is not required and that the keyboard/display would be replaced by 
another interface component, as needed by the application. 

Please replace the paragraph at page 47, lines 23-27 with the following rewritten paragraph: 
We demonstrate the preferred floating ambiguity resolution process using forms [12] - 
[17] above. It may be appreciated that the other previously-described embodiments of the 
present invention may [[use]] also use this floating ambiguity resolution process by simply 



W02-SF:5HY\61424761.1 



-20- 



Application Serial No. 10/670,116 

Amendment Dated August 16, 2004 

Reply to Notice of Allowance of May 14, 2004 

omitting forms and/or terms of forms from the process. The general view of the preferred 
ambiguity resolution process is to reduce the value of the following form during a series of 
epochs: 

Please replace the paragraph at page 48, lines 6-10 with the following rewritten paragraph: 
In form F% , oflfy represent represents the combined vector of current estimates of floating 

ambiguities: A 1>0 N^ = [ A 10 N L1 , A 1?0 N L2 ], and the following are used for the cases of q=l : 
A q , 0 bI L1 = A q , 0 b^, A q , 0 b; L2 =A q?0 b^ 2 , A q50P ; L1 = A q , 0 P^ 5 and A^ 2 =A q , oP ^ 2 . We 
have also broken down the weighting matrices [c„ 1 ]~ l , [cy 1 j" 1 , etc., as follows: 

Please replace the paragraph starting at page 49, line 26 and ending at page 50, line 10 with 
the following rewritten paragraph: 

An exemplary estimation process for the floating ambiguities according to the present 
invention employs form [31] in an iterative manner. We start at the initial epoch k=0 with the 
weighting matrix Do set to the zero matrix, and an initial guess of floating ambiguities A^qNq 
equal to zero. The first term of form [31] evaluates to zero for this initial epoch. We then 
generate a set of values for 6X k , k\$?k > 8 1,0 I A: > A 1,0 N * at a first epoch k=\ which moves the 
value of F s towards zero. This generates an initial estimate A l5 o N i f° r the floating ambiguities 
and the rover's position (by way of 6Xi ). A weighting matrix Di for the ambiguities is then 
generated, and a new set of values is generated for 8X^ , A 10 ^ > ^lja^k > ^Ifl^k at the next 
epoch k=2 which moves the value of F 2 (k=2) towards zero. As a result, subsequent estimates 
^1,0^2 SX 2 are generated. This iteration process continues, with Ax^N* and 8X^ 
generally improving in accuracy as the iterations progress. 
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Please replace the paragraph at page 55, lines 6-8 with the following rewritten paragraph: 
The means for performing [[then]] the above steps in rover 100 are provided by main 
processor 110 under the direction of instruction sets stored in memory 1 12, with the various 
computed data being stored in data memory 114. 

Please replace the paragraph starting at page 56, line 5 and ending at page 57, line 1 with the 
following rewritten paragraph: 

The columns are modified by substituting a value of -1 for each column element that would 
normally be zero in the identity matrix. The positions of columns of matrix E which are 
associated with satellites S x and S 2 correspond to the row positions of these satellites in the 
vector Aj qN^ . When matrix E is multiplied onto vector A^qN^ , the floating ambiguities 

corresponding to satellites s and s remain unchanged, but the floating ambiguity associated 

with s Ll is subtracted from the other floating ambiguities in the LI -band, and the floating 

ambiguity associated with s L1 is subtracted from the other floating ambiguities in the L2-band. 
Next, a permutation matrix II is generated, and is applied to matrix E to gen e rated generate a 
matrix product II -E. Permutation The permutation matrix is constructed to move the S ] and 
S 2 columns of matrix E to the first and second column positions of the matrix product II -E. The 
construction of permutation matrices is well known to the field of mathematics. A permutation 
matrix satisfies the following relationships: II T n = IMI T = I. Next, a change of variables for 
the cost function F is undertaken as follows: 

Please replace the paragraph at page 58, line 17-19 with the following rewritten paragraph: 
Having resolved the fixed-integer ambiguities, we can [[new]] how generate a further 
refined estimate of the other solvable unknowns: 8X^ , A 10 r^ , S^I* • A third cost function is 

formed as follows: . 
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Please replace the paragraph at page 59, line 8-10 with the following rewritten paragraph: 
The means for performing [[then]] the above steps in rover 1 10 are provided by main 
processor 110 under the direction of instruction sets stored in memory 112, with the various 
computed data being stored in data memory 114. 

Please replace the paragraph at page 60, lines 7-13 with the following rewritten paragraph: 
With this information, the residuals (difference quantities) L q r b^ lyS , /S. q ^ r b^ 2 ' s , 

A q r p^ s , and A q r p^ 2,s of forms [4A-4D] are generated for each baseline (q, r), where (q, r) 
has the following pairings (B2, Bl), (B3, Bl), and (B2, B3). For each baseline, the solvable 
unknowns in forms [5A-5D] are A^ r r£ , A qr l£ > m d ^q.r^k • Values are estimated from the 

residuals in a manner similar to that described [[about]] above for the primary baseline between 
the Rover and first base station. As an example, a cost function F(*) may be formed as follows: 

Please replace the paragraph at page 61, lines 3-14 with the following rewritten paragraph: 

The floating ambiguities may be estimated in an iterative manner as described above. 
We start at the initial epoch £=0 with the weighting matrix D 0 set to the zero matrix, and an 
initial guess of floating ambiguities A q?r N 0 equal to zero. The first term of form [31] evaluates 

to zero for this initial epoch. We then generate a set of values for A^ r r^ , Aq 5r I^ , Aq 9l . at a 
first epoch k=l which moves the value of F towards zero. This generates an initial estimate 
A q,r Nj for the floating ambiguities. A weighting matrix Di for the ambiguities is then 
generated, and a new set of values is generated for kq^k > Aq,rljt > ^q^k at the next epoch 
k=2 which moves the value of F (k=2) towards zero. As a result, a subsequent estimate s 
estimate A q5r N 2 is generated. This iteration process continues, with A^N^ generally 

improving in accuracy as the iterations progress. The corresponding H*Y = B form for this 
process (similar to form [32]) is as follows: 



W02-SF:5HY\61424761.1 



-23- 



Application Serial No. 10/670,116 

Amendment Dated August 16, 2004 

Reply to Notice of Allowance of May 14, 2004 

Please replace the paragraph at page 64, lines 6-9 with the following rewritten paragraph: 

The corresponding fixed-integer ambiguities A q?r N^ may be generated from the floating 

ambiguities Aq^N^ by the same process described above with reference to forms [37] - [43] 

that [[is]] was used to generate the fixed-integer ambiguities associated with the primary base 
Ime baseline between the rover (R) and first base station (Bl). 

Please replace the paragraph at page 64, lines 10-12 with the following rewritten paragraph: 

Having resolved the fixed-integer ambiguities, we can [[new]] now generate a further 
refined estimate of the other solvable unknowns: A^ r r^ and A q>r l£ . A third cost function is 

formed as follows: 

Please replace the paragraph starting at page 64, line 10 and ending at page 65, line 8 with the 
% following rewritten paragraph: 

where the five terms of the form are the same as the second through sixth terms of form [45], 
except that floating ambiguities have been replaced by the fixed-integer ambiguities. The same 
estimation procedures used on form [45] may be applied above to form [51], except that only one 
iteration is needed, and matrix D is not generated or used. As a result, estimates for A q ^ r T k and 

Aq 5r I^ are generated. However, we prefer to perform some consistency checks on these 

estimates before providing them to the process that operates on the primary baseline between the 
rover and the first base station. Thus, we will denote these estimates with hat symbols as 

follows: A^ r f £ and A q r i^ . After the estimates associated with the three ba se line baselines 

(q,r) = (B2, Bl), (B3, Bl), (B3, Bl) have been generated for the k-th epoch, we have the 
following data: 
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Please replace the paragraph at page 65, lines 17-23 with the following rewritten paragraph: 

If the above relationships are not satisfied, the ambiguities have been resolved incorrectly for at 
least one ba s e lin e baselines , and the estimations of fixed ambiguity should be neglected. To 
address an incorrect resolution, new data may be taken, or various subsets of data may be 
proc e s s e s processed to generate sets of ambiguities which satisfy the above relationships. After 
all three between ba s e rec e iver s ambiguity between-base receivers 1 ambiguities have been 

resolved and the above relationships are satisfied, the ambiguity vectors &2,l^k > ^3,1^ , and 
^3,2 ar e considered to be correctly fixed. 

Please replace the paragraph at page 66, lines 4-12 with the following rewritten paragraph: 

Next, we perform a consistency check on the ionosphere delays. In the following 
discussion, [[we]] the subscripts on the ionosphere delay differentials have sometimes been 

exchanged, but this is of no substantive consequence since A 1?2 I# = -^24!* 5 

^2,3!* = -^3j\k> ^1,3! k = ~ A 3 5l l£. The between base rec e iv e rs between-base 
receivers 1 ionosphere estimations should satisfy the relationship 

A 12 I ic + k + A34I ^ = 0 , to within a tolerance value of ±e 2 . [54] 

However, mea s urement s measurement noise typically prevents these relationships from being 
satisfied to an acceptable tolerance level ±6 2 . To better satisfy [54], the following quadratic 
function is minimized to obtain new estimations A 12 1* > ^3 Ifc , A34 1# : 

Please replace the paragraph starting at page 66, line 22 and ending at page 67, line 1 with the 
following rewritten paragraph: 

Before applying form [57], the vectors A^I* > A 2,3^ > ^3 9 i^k ^ preferably smoothed using a 
Kalman filtering scheme with a dynamic model based on the above de s cribed above-described 
[[the]] Gauss-Markov time model (forms [GM1] and [GM2] above). 
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Please replace the paragraph at page 67, lines 2-8 with the following rewritten paragraph: 

The results of these proce ss processes may be provided to the process which generates 
the ionosphere delay differentials for the baselines associated with the rover stations, specifically 
the differentials generated according to forms [26], [29], and [30]. For form [26], which 

generates A^qI^ , the results are provided as follows: A^I^ = A^I* ^1,3^ k = Ai^I* • 
For generating the ionosphere delay differentials according to form [29] w e [291, we set 

= ~^\,2^k ■ For generating the ionosphere delay differentials according to form [30] we 
[30], we set A 34 I* = A 3>1 T it . 

Please replace the paragraph starting at page 67, line 29 and ending at page 68, line 3 with the 
following rewritten paragraph: 

While it is preferable to interpolate the ionosphere delays using two baselines between 
three base stations, it may be appreciated that some application applications may achieve 
acceptable accuracy by [[only]] interpolating the ionosphere delays using only one baseline 
between two base stations. Such an example may be a road project where the road is relatively 
straight, as shown in FIG. 4. 

Please replace the paragraph at page 68, lines 4-9 with the following rewritten paragraph: 

Each of the above methods of generating the base station data and estimating the 
coordinates of the rover is preferably implemented by a data processing system, such as a 
microcomputer, operating under the direction of a set of instructions stored in computer-readable 
medium, such as ROM, RAM, magnetic tape, magnetic disk, etc. All the methods may be 
implement s implemented on one data processor, or they may be divided among two or more 
data processors. 
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Please replace the Abstract of the Disclosure at page 100, lines 5-16 with the following 
rewritten Abstract: 

Disclosed are methods and apparatuses for determining the position of a roving receiver 
in a coordinate system using at least two base-station receivers, which are located at fixed and 
known positions within the coordinate system. The knowledge of the precise locations of the 
base-station receivers makes it possible to better account for one or all of carrier ambiguities, 
receiver time offsets, and atmospheric effects encountered by the rover receiver, and to thereby 
increase the accuracy of the estimated receiver position of the rover. Base line s Baselines are 
established between the rover and each base-station, and a ba s e line is baselines are established 
between the base stations. Navigation equations, which have known quantities and unknown 
quantities, are established for each baseline. Unknowns for the baseline between base stations 
are estimated, and then used to correlate and reduce the number of unknowns associated with 
rover baselines, thereby improving accuracy of the rover's estimated position. 
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